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Abstract: Integration By Parts (IBP) is an important method for computing Feynman 
integrals. This work describes a formulation of the theory involving a set of differential 
equations in parameter space, and especially the definition and study of an associated 
Lie group G. The group acts on parameter space and foliates it into G-orbits. The 
differential equations essentially provide the gradient of the integral within the orbit in 
terms of integrals associated with degenerate diagrams. In this way the computation of 
a Feynman integral at a general point in parameter space is reduced to the evaluation of 
the Feynman integral at some freely chosen base point on the same orbit, together with 
a line integral inside the G-orbit and the degenerate integrals along the path. This paper 
restricts to vacuum diagrams and integrals without numerators, but the method is expected 
to generalize. The method is demonstrated by application to the two-loop vacuum diagram. 
A relation between the reducibility of a diagram though IBP and the reducibility of the 
associated electric circuit though the Y — A transform is speculated. 
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To the reader 


1 Introduction 

Integration By Parts (IBP) is an important method for computing Feynman integrals, 
discovered in 1981 [1]. It should be distinguished from the general method of effecting 
integration which bears the same name. IBP enables to compute many diagrams, for 
instance the 2-loop propagator diagram shown in fig.l, and whenever it applies it works 
for all space-time dimensions. Ref. [1] currently has about one thousand InSpire citations, 
evidence for both its wide application and for its theoretical interest. 

Despite its popularity and 34 years of research, so far it is not known to which diagrams 
the method should be applied and what output should be expected of it. Already the 
original authors remarked (in the discussion section) “As yet there is no criterion to decide 



Figure 1. The two-loop propagator diagram. The associated Feynman integral is conveniently 
computed by the IBP method. 
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whether a given 4 loop p-integral [propagator type integral - BK] can reduced to simpler 
ones via integration by part” as well as “All these observations leave the impression of 
something important having been missed in our analysis, that could be very useful for both 
practical purposes and a better understanding of perturbative series, if these two things 
can be separated” . Indeed this question is natural for any method, a sort of ‘declaration 
of purpose’ to appear in the ‘instruction manual’. (It is hoped that the quotations would 
not trigger the arXiv ‘text overlap’ alert...) 

We can mention here a few more specific indicators for the incompleteness of the current 
theory for the IBP method. Presently the practice is to use Laporta’s algorithm [2] realized 
through computer packages both for obtaining the IBP recursion relations, and then for 
solving the recursion. Indeed, computers are very useful to carry long computations, yet 
blind usage of ‘black boxes’ may obscure the general properties of the method. In addition, 
the very name of the method, I find to be lacking, since it shares the name of the general 
method of integration, while the potential of IBP depends strongly on the diagram topology, 
and that is not reflected by the name. This is another indication for an incomplete state 
of understanding. 

In this paper I shall carry some steps to discover the ‘statement of purpose’, and I 
believe the full answer will be found soon. Section 2 defines a generating function for Feyn¬ 
man integrals for a given diagram topology, and its computation would be our goal. Section 
3 describes how to obtain IBP relations both in index space n and in mass-squared space 
X. The main novelty appears in section 4 where we motivate and define the numerator-free 
sub-algebra of IBP relations, which leads to a system of linear first order partial differential 
equations. We proceed to outline the general procedure of solving the equations in section 
5. Next in section 6 we demonstrate the ideas of this paper through a concrete diagram - 
the two loop vacuum diagram. Finally we offer a summary and discussion in section 7. 

Being a newcomer to the subject the author may have omitted relevant references, and 
he would be happy to receive correspondence on this and other matters. 

Background. The IBP method was discovered in [1]. Baikov’s method, a (non-recursive) 
solution to the recursion relations, was presented in 1996 [3]. During the 90’s the closely 
related method of Differential Equations (DE) was developed [4-6], see the reviews [7, 8] 
where the last one also reviews recent progress. The IBP and DE methods occupy a 
prominent place in the excellent textbooks [9, 10] on Eeynman integrals. Fairly recent works 
on IBP include: an introduction to IBP [11]; a comparison of master integral counting 
between IBP and DE [12]; an algorithm for the generation of unitarity-compatible IBP 
Relations [13]; new relations between master integrals using IBP [14]; a recent version of 
IBP computerized tools [15]; an improvement over Laporta’s algorithm which employs a 
Monte-Carlo approach [16]; a suggested improvement for Laporta’s algorithm using finite 
fields [17]; a recent work from the viewpoint of differential geometry utilizing differential 
forms [18]; and finally an application to multiploop Four Dimensional Regularized integrals 
[19]. 

The author learned about the IBP method through its application to the analysis of 
the post-Newtonian limit of the two-body problem in Einstein’s gravity [20, 21] (which is 
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a non-quantum problem). 


2 Generating fnnction 

Given a Feynman diagram F we would like to compute a general integral of the form 

where dl ;= 0^=1 denotes integration over all loops for a general space-time dimension 
d; the index i = 1,... ,P runs over all propagators; is the energy-momentum of the 
propagator and it should be considered to be expressed as a combination of loop and 
external momenta; the x parameters are standard notation and can be identified with 
squared masses Xi = m?; and finally the powers n* are termed indices. /r(a^,n) is defined 
as a very general Feynman integral associated with F. The indices Uj generalize the common 
case n* = 1 and are necessary as they may appear during IBP manipulations of unit indices. 
More generally the integral could also include numerators with arbitrary powers of Lorentz 
scalars, but for the current purposes (2.1) suffices. 

The general integrals Ir{x, n) for all natural indices n* = 1, 2, 3,... can be conveniently 
encoded into the generating function 

Ir{x) := I [kl-x ) ^ ^r(a:,n)|„,=i (2.2) 

Indeed I{x,n) can be recovered through 

Since I{x,n) is encoded in the Taylor coefficient of I{x) it is appropriate to refer to the 
latter as a generating function. Note that whereas an ordinary generating function encodes 
a single series of coefficients through its Taylor coefficients at one specific point, such as the 
origin, here the Taylor coefficients of I{x) around any point are of interest. In addition, 
since I{x) = I{x,ni = 1) this generating function has a very clear interpretation, and the 
x variables carry a double interpretation - both as formal parameters of the generating 
function as well as mass squares. Hence we find that the integrals I{x) contain as much 
information as I{x,n) and we set a goal to compute /r(a^)- 

3 IBP relations and generators 

The standard way to derive the IBP relations is to consider the identity 

0 = jdl^^kl (3.1) 
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where d/dl = d/dl^ is a divergence with respect to one of the loop momenta; k = 
is a propagator momentum and I is any integrand, which we will mostly take to be the 
integrand of (2.1), namely 

Ui {kf - Xi) 

The vector indices n are suppressed for clarity, and more generally k could be a linear 
combination of propagator momenta, or equivalently of loop and external momenta. 

By Gauss’s divergence theorem the RHS of (3.1) equals a flux integral over a surface at 
infinity. For small enough dimensions this integral vanishes (in the limit that the surface 
approaches infinity) and hence in dimensional regularization, which is assumed here, it 
vanishes identically. The method’s name is derived from the use of the divergence theorem 
which is a generalization of the elementary integral identity of integration by parts. [1] adds 
another point of view, that of configuration space, explaining that the divergence identity 
(3.1) is related to the vanishing sum of 3 vectors representing triangle edges (eq. 3.4b over 
there). 

In order to obtain the IBP relations we proceed to expand the RHS of (3.1). When 
k lies in the loop I its differentiation generates a term with a factor of d. Next we should 
differentiate the propagators in the denominator. For any ki lying in the loop I we get a 
term of the form 

2k ■ ki rii i"'' (3-3) 

where 1"*“ is a raising operator which increases the integrand’s index m by 1. 

Now comes a crucial point. The Lorentz scalar k ■ ki can be expanded as 

2k-h = Y,Tik‘] + M, , (3.4) 

j 

namely into a combination of squares kj plus possibly a remainder, denoted M*. Such 
factors can appear in the numerators of more general Feynman integrals. For example, in 
a trivalent junction where k 2 = k — ki we have 

2k- ki = + k\ — k2 (3.5) 

Defining a propagator term 

E := k'^ — X (3.6) 

we can replace /c| = Ej + Xj and then we further replace Ej —>• j~, where j“ is a lowering 
operator which reduces rij by 1. Finally we arrive at the following types of terms which 
can appear upon expanding the RHS of (3.1) 

d , nii^{]~+Xj) , nji+Mj (3.7) 

The operators i^, i“ relate integrals with different values of n* and hence we obtained 
a recursion relation for I. This relation is clearly linear in /, and it is of order at most 1 
both in raising and in lowering operators. In fact, by choosing different values for I and k 
we obtain a set of linear recursion relations, known as the IBP relations. More specifically 
for a vacuum diagram with L loops we may generating an independent set of equations by 
choosing I = k, k = Ig where r, s = 1,..., L run over all loops. 
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3.1 IBP relations in x space 


We notice that the terms appearing in the IBP recursion relations (3.7) are special. A 
term of the form /(nj) i^j~ with /(nj) an arbitrary function of n* would preserve linearity 
and the order of the IBP recursion relation. However, we notice that the dependence on n 
appears in a special form n* i"*", namely, a factor of is always accompanied by a raising 
operator. 

This leads us to seek a set of variables where this property may have a more natural 
interpretation. In general it is known that recursion relations are connected with differential 
equations: a differential equation becomes a recursion equation once we substitute in its 
Taylor series expansion. Conversely, a recursion relation can be equivalent to a differential 
equation for its generating function, in some formal parameters x. In this way we could 
motivate the definition of the generating function (2.2) even if we started with completely 
massless integrals in n space, namely I{xi = 0,ni). In fact one can consider various ways 
of collecting a series I{n) into a generating function I{x) := f i’’^) where 

/(n) could be an arbitrary function of n, such as / = 1 or / = l/(re — 1)!. The choice (2.2) 
corresponds to choosing /(n) = 1 and is exactly such that the factor rii i"'' is replaced by a 
derivative dfdxi. 

An IBP relation in x space is derived by starting with the same identity (3.1) where 
this time all indices are set to unity in the integrand I (3.2). Whereas in the recursion 
relation context the mass squares Xi were considered fixed, here we notice that we can 
replace riii^ —)• djdxi as anticipated. The appearance of djdxi can be thought to arise 
from the identity 

if = -i; 

which holds for all functions of the form I = I{k‘j — Xi). Now the derivative djdxi can 
be taken outside the integral, and then an IBP relation can be stated at the level of the 
integral I{x) as follows 

0 = [c + T) Xi d^] I + J (3.9) 

where c and T^ are constants (c may depend on d as follows c{d) = cid — C 2 ) , = d/dxj 

and finally J is a source term 

J = j dl [Tj Ei + Mj cP] I 
where Mj are some numerators. Noticing that 


j dlEil = Oil 


(3.10) 


where Oj/r := lOiV and OjP is the diagram obtained from P by omitting propagator i, the 
first term of J can be written in terms of generating functions for degenerated diagrams, 
namely 


J = Tj Oil + d^ dl Mj I . 


(3.11) 
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The differential equations appear to be the same at those suggested by [22], and are 
also related to the approach of [3]. 

Altogether in x space we arrive at the differential equation (3.9) with the source term 
(3.11). If there are several IBP relations labelled by a then we obtain a set of differential 
equations and their parameters obtain an a index as follows 


0 = 


c“ + (r“)'. 


/ + 


(3.12) 


where depends now on M“. 

These equations are linear in I and first order in the derivatives djdxi. It should be 
stressed that the first degree in derivatives is not a direct consequence of the first degree 
of the recursion relations, but rather it relies on the special dependence on n* which we 
noticed. In this sense the IBP relations are more natural in the x plane. 


4 Lie algebra structure and sub-algebras 

The IBP equations in x-space (3.12) immediately suggest a Lie algebra structure (it is 
enough that the equations are linear and of first order in derivatives). It is convenient to 
define differential operators 

:= c“ + (r“)V X, (4.1) 

such that the differential equations become (3.12) 


0 = I + 


(4.2) 


Consider the commutator [D°‘,D^]. The terms quadratic in derivatives cancel, so the 
commutator is linear in derivatives. As such it generates another equation of the same 
type. If we already had all the equations, then this commutator must be linear combination 
of the differential operators, namely 




= . 


(4.3) 


This is a Lie algebra structure with structure constants 

We can be a bit more explicit about the commutation relations of the differential 
operators (4.1) 

Xi (4.4) 

j 

where the commutator on the LHS is the operator commutator and on the RHS it is simply 
the matrix commutator. 

The Lie algebra structure of the IBP relations was described already by R. Lee in [23] . 
There the algebraic structure was studied in relation to the Lorentz Invariance identities 
for external momenta. It was described in n space and studied through the commutators 
of the IBP generators Ors = d/dlrh^ which are the same operators appearing in (3.1). 
Clearly the group structure among the differential operators (4.3) is inherited from 
the commutation relations of the associated operators (namely, the mapping 0—7-44 
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preserves the Lie algebra structure). That paper does not mention the x space form of the 
IBP equations, where as we have shown the Lie group structure is apparent. 

Applying the commutator to the IBP equations (4.2) we obtain a constraint. We start 
with 


0 = [D^I + J^] -D^ (D“I +J“) = 




I + -D^ / 


(4.5) 


Using (4.3) in the last line and comparing with 0 = I + J'^) we now obtain 


jb _jja jb^ jab^j, 


a jb 


fab TC 


(4.6) 


This is an integrability constraint which the sources must satisfy. 


4.1 The numerator-free sub-algebra 

Defintion. We denote by Gq the Lie algebra of all IBP relations, namely all the operators 
(4.1) (and the associated equations (3.12) ) obtained from all possible identities of the 
form (3.1) for arbitrary I and k vectors (in the presence of external legs these identities 
should be supplemented [3, 24]). 

Considering the source term (3.11) we observe that the first term contains degenera¬ 
tions of P and hence will be assumed to be known (otherwise we have no right to expect to 
be able to compute Ir) while the second term includes a more complicated and therefore 
possibly unknown integral. This motivates the following definition: 


Defiuitiou. We denote by G the subspace of Gq such that all sources (3.11) are free of 
numerators. 

A similar motivation can be seen in n space. In the presence of numerators the IBP 
relations contains terms which increase an index in the numerator (as well as the index of 
a denominator), and hence the resulting recursion relation could be considered ineffective, 
as the integral with the numerator would be usually unknown. 

G is not only a sub-space but rather a sub-algebra of Go, as the commutator of two 
numerator-free relations is again numerator-free. We refer to G as the numerator-free Lie 
algebra. 

From now on we restrict ourselves to vacuum diagrams and integrals without numer¬ 
ators (2.2). Including numerators is possible, but as we have seen it is necessary first to 
solve for the numerator-free integral. Including diagrams with external legs complicates 
the discussion a bit, and so incorporating it in the current formulation is reserved for future 
work, though in principle it is understood [3, 24]). 

Let us summarize the set of equations for I(x) associated with G 

0 = D“I-kJ“ (4.7) 

where 

T)“ = c“ + {T^)G Xi 

ja ^ gj qj 
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where a runs over the generators of G. This is a set of linear hrst order PDE’s. Moreover, 
the operators (4.1) are homogeneous in the x variables. 

The definition and study of G and its associated equation set are one of the central 
results of this work. The suggested approach is to associate with any Feynman diagram T its 
numerator-free Lie Algebra, to study it and use it to appropriately reduce the corresponding 
Feynman integral as much as possible. 

The name Integration By Parts is misleading from the perspective of G since integration 
by parts is a general technique which applies to any integral whereas here the requirement 
to avoid numerators is essential and it reflect the diagram’s topology. 

We note that [25] showed that numerators can be eliminated by raising the dimension. 

Constant-free sub-algebra 

Another useful sub-algebra is given by the following definition: 

Definition. We denote by G 2 the subspace of G such which is constant-free, namely such 
that c“ = 0 in (4.7). 

G 2 is co-dimension 1 in G, it is a sub-algebra of G and in fact an ideal, namely 
[ 9 ^ 92 ] £ G 2 for all g G G, g 2 ^ G 2 - For a vacuum diagram without numerators, the 
complement to G 2 is the equation 

{) = DQl=[cQ-Xid^]l (4.8) 

where cq is the x dimension of I. In other words, this equation is the Fuler theorem for 
homogeneous functions reflecting dimensional analysis of I. We call Dq := cq — Xj 9* the 
dimension differential operator. 

5 Outline for method of solution 

Given a vacuum Feynman diagram P and a numerator-free integral (2.2) we associated 
with it a set of equations (4.7) and a numerator-free Lie algebra of differential operators 
(4.3). In this section we outline the method of solution for this set. We reiterate that the 
current formulation is expected to generalize further to allow for external legs and integrals 
with numerators. 

Any single equation or group generator, namely a specific index a in the equation set, 
defines the standard field of characteristic curves Xj = Xj{s) through 

^Xj = {T^Yj Xi (5.1) 

These are first order equations, which in this case are simply linear with constant coeffi¬ 
cients. 

In the presence of a set of equations we have a set of characteristic curve helds. This 
set generates characteristic sub-manifolds which foliate the x space (the co-dimension of 
the sub-manifolds could be greater than zero as a result of the closure of the commutation 
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relations). In fact these sub-manifolds are nothing but orbits of the symmetry group G in 
X space. 

It is useful to imagine coordinates which are adapted to the group action: G orbits are 
parameterized by while Pi,P 2 , ■ ■ ■ are independent invariants of G and can be used as 
coordinates for the transverse direction. 

The equation set (4.7) can be solved algebraically to yield the gradient of I within the 
G-orbits as a function of I and x: d/dE'^I = /“(/, x). Therefore the computation of I{x) 
can be reduced as follows 


• Evaluate I = I{Pi,P 2 ,...) on a freely chosen sub-manifold Xq of initial conditions 
which is transversal to the G-orbits. 


• To evaluate I{x) for any x, hnd xq G Xq on the same G-orbit as x, connect them by 
a curve x*^ = x^(s), such that x*^(0) = xq, x^(1) = x and hnally solve an ordinary 
differential equation along the curve, namely 


A 

ds 


i = r{i,x) 


ds 


(5.2) 


Clearly the co-dimension of a G-orbit is central to this approach as it sets the dimension 
of the initial value surface Xq. The freedom in choosing Xq can be exploited to choose 
integrals with a higher degree of symmetry, for instance, several identical masses (x’s). 

We can reduce the problem a bit further. Let us solve the homogeneous version of 

(4.7) , namely 

0 = lo (5.3) 

If we consider only generators lying within the constant-free subgroup G 2 then by definition 

Iq is constant on G 2 orbits. This leads us to split the coordinates into coordinates |q, 
which parameterize the G 2 orbit and one other coordinate, which we denote by P. Then 
we should consider the the dimension equation Dq Iq = 0. It is an ODE in the dependent 
variable P, implying some solution Iq = c(Pi, P 2 ,...) Io{P)- 

Now solutions for the original, inhomogenous equations can be obtained as usual 
through the method of variation of the constants 

/(x) = c(x) Jo . (5.4) 

Upon substitution back into (4.7) one hnds the projected gradient of c, d/df^aC = /“(x) 
where as usual /" is now independent of c and hence the solution for c reduces to a line 
integral 

c{x) = c{xq)-\- [ p{x)dia (5.5) 

Jio 

Finally I would like to mention another possible approach to solving system of equations 

(4.7) over the G-orbits. One can consider not only first order differential equations, but 
also higher order differential equations. For example, one could consider a second order 
PDE formed from a quadratic Casimir of the differential operators D°^. 
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Figure 2. Trivalent vacuum diagrams with few loops, (a) The 2-loop vacuum diagram. Its 
associated Feynman diagram is analyzed in the text, (b) The 3-loop trivalent vacuum diagram. It 
has a tetrahedron topology 


6 Demonstration 


In this section we shall analyze a specific Feynman integral with the method described 
above. But hrst we start with a more general perspective. In this paper we consider 
mostly vacuum diagrams and associated integral without numerators. It is illuminating to 
find the relation between the number of numerator terms and the number of loops. 


Number of numerators 


A vacuum diagram with L loops defines Q = L(L-|-l)/2 independent quadratic expressions 
(Lorentz scalars). The number of propagators P, on the other hand, assuming trivalent 
vertices is P = 3(L — 1) [this is gotten from the following two equations: the Euler char¬ 
acteristic gives L — P + V = 1 where V is the number of vertices and the trivalent nature 
of the vertices implies 3V = 2P]. Hence the number of numerator terms is 

#Num = Q-P= _ 3(L - 1) = 1(L - 2)(L - 3) (6.1) 


Let us tabulate this function for some low values of L 


L 

2 

3 

4 

5 

#Num 

0 

0 

1 

3 


The IBP relations of a vacuum diagram are generated by the operators 


which generate the algebra 


Ors = d/dlr Is 


Go = GL{L,R) 


of invertible real L by L matrices. 

We shall now discuss some cases one by one. 


(6.3) 

(6.4) 
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There is a single 2-loop trivalent vacuum diagram, see fig. 2(a). There are no numer¬ 
ators in this case, so all IBP relations are numerator-free and 

G = Go = GL{2,R) . (6.5) 

At 3-loops there is still a single such diagram ~ the tetrahedron, see fig. 2(b). In this 
case again all IBP relations are numerator-free and so 

G = Go = GL{3,R) . (6.6) 

At 4-loops there are already 2 inequivalent trivalent vacuum diagrams, and there is a 
possible numerator. So we can expect that G would be a proper subgroup of GL(4,M). 

Each non-vacuum diagram defines its “vacuum closure” in the usual way - by attaching 
all external legs to a new point at infinity. 


Two-loop vacuum diagram 

We now turn to demonstrate the method by applying it to a specific diagram, the 2-loop 
(trivalent) vacuum diagram of fig. 2(a). It was computed in [26] using the Mellin-Barnes 
transform, see also references therein. 

The symmetry group of the diagram allows to exchange the two vertices and to permute 
the 3 propagators and hence has 12 elements. 

There are 3 propagators, so the generating function (2.2) depends on 3 x variables 
(mass-squares) xi, X 2 , X 3 . According to table 6.2 there are no possible numerators, and 
so this function generates all the possible Feynman integrals. For a specific choice of loop 
currents its definition becomes 

^ f d%d% 

I [Xi, X2, X3) I , 2 \ /,2 \ //7 I 7 ^2 ^ 

j [li - xi) - X2) {{h + hr - X3) 

The set of all possible IBP’s generates the group Go = GL{2,R) according to (6.4) 
with L = 2. As there are no possible numerators, there are no possible obstructions for 
the definition of the numerator-free sub-algebra, and we have G = Go = GL{ 2 ,R). 

In order to form the equation set (4.7) we need to know the Feynman integrals of all 
possible degenerations. In this case there is a single degeneration into the figure 8 diagram 
shown in fig. 3, whose value factorizes 

Icx,{xi,X 2 ) = jixi) j{x 2 ) (6.8) 


and can be considered to define the function j. 

The equation set can now be cast in the following rather symmetric form 


d — 3 —xi —X2 —X3 


' 1 ' 


0 

0 X2 - X3 X2 -X3 


di 

1 + 

ih - 22 ) fi 

0 —xi X3 — xi X3 


82 

Ui - js) 32 

0 Xi — X2 Xi — X 2 


1 

_ 1 


_ {h - h) 3 3 _ 


(6.9) 
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Figure 3. The one and only possible degeneration of fig. 2(a). 


where the functions ji which appear in the sources are given by 

ji-=j{xi) ( 6 . 10 ) 

The first row can be identified with the dimension generator Dq (4.8), while the next 
three rows generate the group G 2 = 5'L(2,M) C G. More explicitly, the first generator is 

Li = {x2 - X3)d^ + X2d'^ - ( 6 - 11 ) 


and ^ 2,^3 are defined by cyclic permutations. The commutation relations among the Lj’s 
are 


[Li,L 2] — —Li — L 2 


( 6 . 12 ) 


and the rest can be obtained by cyclic permutations. 

We confirmed that the sources satisfy the integrability constraint (4.6). 

The X parameters transform in the same way as the quadratics Ir'h: so they are in the 
symmetric representation, which indeed has dimension 3. In fact we can define a symmetric 
matrix in terms of x 


X'ps 


2xi -Xl - X 2 + X^ 

-Xi - X2 + x^ 2X2 


(6.13) 


then g £ G = GL{2, M) acts on it through 


Xrs 


^ 9 Xrs 9 


T 


(6.14) 


The orbits of G are 3-dimensional so the initial conditions sub-manifold is 0-dimensional, 
namely a point. In fact, the integral is known at the following two points xi = X 2 = 0, X 3 = 
w? (e.g. fig 3.4 on p.37 and eq. (10.39) of [10]) as well as xi = 0, X 2 = X 3 = w? (e.g. 
example 3.4 on p. 48 and eq. (10.38) of [10]). 

In order to obtain homogenous solutions for (6.9) we notice that there is a single G 2 
invariant x expression, namely 

P := det {xrs) = {xi + X 2 + ^ 3 )^ -2{x\ + X 2 + x|) (6.15) 

We notice that P is symmetric in Xi, as it should, even though the matrix Xrs was not 
so, due to the breaking of the symmetry by the choice of loop currents. This is the same 
polynomial that appears in Baikov’s method [3] when applied to this case, see [9] eq. (6.17) 
after the replacement —)■ X 3 . 
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( 6 . 16 ) 


The first row of (6.9) allows now to solve for Iq to obtain 

Jo = const 

Defining 

I{x) = c{x)Io{P) (6.17) 

we can solve (6.9) to obtain the gradient of c and thereby reduce the solution of I{x) to a 
line integral. 

Alternatively we may consider the 2nd order differential equation given by the quadratic 
Casimir 

A = (Li + L 2 + T 3 ) — (T^ + + L|) (6.18) 

7 Summary and Discussion 

Summary. Our goal is to fully exploit the IBP method to compute Feynman integrals. 
The standard method is to employ recursion relations to reduce the given integral to a 
set of master integrals using a computerized algorithm and then to compute the master 
integrals, possibly through the method of differential equations. 

We find that the generating function I = I{x) defined in (2.2) such that x denotes mass- 
squared, satisfies a set of differential equations (4.7) in x space indexed by the generators of 
a Lie group G. The group G appears to be a new concept. It depends only the diagram’s 
topology G = Gr, where T represents the diagram, and can be used to characterize it. 
The group G foliates x space into G orbits. The PDF set reduces the computation of the 
general Feynman integral to the evaluation of the Feynman integral over an initial value 
sub-manifold of parameters which is transversal to the G-orbits, followed by the evaluation 
of a line integral. Another approach using higher order differential equations was suggested. 

The method was demonstrated for the case of the 2-loop vacuum diagram. This 
diagram consists of 3 propagators and hence we have 3 x’s: xi, X 2 , X 3 . The initial group 
Go is Go = GL(L,M) = GL(2,M). All IBP relations are numerator-free so the group 
G = Go = GL( 2 ,]R) and the constant-free group is G 2 = S'L( 2 ,M). The orbits of G are 3d 
(co-dimension 0). Hence the surface of initial conditions is a single point xq to be chosen 
at will. In fact there two are points xq where I{xq) is known. The orbits of G 2 are given 
by the level-sets of the Baikov polynomial P = (xi -|- X 2 -|- xs)^ — 2(xf + x"^ + ^^ 3 ). The 
homogeneous solutions are Iq = const . In this way the evaluation of /(x) for 

general x is reduced to the evaluation at a single point /( xq ) and a line integral. 

7.1 Discussion 

Interpretation of G. The IBP relations define a Lie algebra, which is particularly evident 
in X space. This algebra generates a group which defines the G-orbits in X space. The 
G-orbits describe the reach of the set of differential equations - they enable to replace the 
required point x by any point xq lying on the same G-orbit. 

The ‘instruction manual’ and the T — A transform. So far we did not address 
the original question, namely when the method should be applied and what the expected 
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benefits should be. We now see that the method is more useful the larger G is, and 
more precisely the smaller is the co-dimension of the G-orbits. There remains the general 
question how G is to be determined given T. 

I would like to speculate further on the diagrams where the method applies. The classic 
IBP relation is the triangle identity. This can be related to a known concept from the held of 
electric circuits. The Kirchhoff equations are known to determine a solution for any circuit. 
However, in certain cases the solution can be obtained through a stepwise reduction, namely 
when components are connected either in parallel or in series. In fact, there is another 
case where stepwise reduction is possible, namely the case of the triangle which can be 
transformed into a Y shaped sub-circuit through the so-called Y — A transform. We would 
like to conjecture that a Feynman integral (without numerators) can be fully solved in n 
space and expressed in terms of T functions alone exactly when the corresponding circuit 
is reducible using the Y — A transform. For example, out of the three 3-loop propagator 
diagrams L, M and N only N, the non-planar one cannot be reduced by IBP to F functions, 
and only its associated electric circuit cannot be stepwise reduced not even through the 
use of the Y — A transform. 

Relation of IBP and DE methods. The methods of IBP and DE and known to be 
related. From the current perspective, they in fact appear to be one, as the IBP relations 
are represented by differential equations in the space of x parameters. 

Relations with the literatnre. In our demonstrated example the homogenous solution 
Io{x) turns out to be the same as Baikov’s function g{x). This is likely not a coincidence. 

In the usual method it is important and interesting to know the number of master 
integrals. In 2010 this number was shown to be finite [27]. In 2013 it was shown in [28] 
to be counted by certain critical points of the Kirchhoff-Symanzik polynomials. Clearly 
this number depends only on the diagram’s topology and one can hope for an even simpler 
method to obtain such a simple integer. This question in fact was a main original motivation 
to the current research. The current results take away emphasis from master integrals and 
their number, and instead suggest to focus on other discrete data, including the identity 
of G and the co-dimension of its G-orbits. 

Open questions. I see many open questions, several of them are currently under study. 
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